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We study the threshold production of two pions induced by neutrinos in nucleon targets. The 
contribution of nucleon, pion and contact terms are calculated using a chiral Lagrangian. The con- 
tribution of the Roper resonance, neglected in earlier studies, has also been taken into account. The 
numerical results for the cross sections are presented and compared with the available experimental 
data. It has been found that in the two pion channels with tv + 7v~ and ir°7v° in the final state, 
the contribution of the TV* (1440) is quite important and could be used to determine the iV*(1440) 
electroweak transition form factors if experimental data with better statistics become available in 
the future. 
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I. INTRODUCTION 



The important ongoing experimental effort addressing the questions of neutrino oscillations is bringing out, as a 
fortunate byproduct, much information on the structure of hadrons and nuclei. Apart from the intrinsic interest of 
the knowledge of axial form factors, structure functions or the strange quark content of the nucleon, a proper and 
precise understanding of various processes induced by neutrino interactions is required in the experimental analysis 
of background substraction, z^-flux determination and particle identification in the neutrino oscillation experiments. 

In the study of the neutrino-nucleus interaction, we can distinguish several energy regions. At low energies we have 
quasiclastic scattering (QE) in which a nucleon is knocked out of the target nucleus. There has been an strong effort 
both experimental and theoretical to measure and describe this process [l|, [2|, y, [J, |5j, [y, [7|, [8|, [9( . The basic ingredient 
is the relatively well known neutrino nucleon elastic interaction, although nuclear effects like Fermi motion, Pauli 
blocking or long range RPA correlations are also needed. These nuclear effects have been found to strongly modify 
the cross section and also to distort angular and energy distributions of the final particles. 

At intermediate energies, above 0.5 GeV, one pion production becomes relevant. The knowled ge o f the elementary 
process v + N — > £. + N' + w is not so well established. One of the reasons is the scarcity of dat a TloL [T l]| . Most of the 
theoretical models assume the dominance of A(1232) resonance mechanisms [HI, H Bill, [II [ISM but others also 
include background terms [jjj The major uncertainties of these models appear in the AA(1232) transition 

axial form factors that are fitted, with some theoretical ansatz to the available data. It is hoped that new data on 
single pion production from neutrino experiments at K2K [12] and MiniBooNE [2^, [24[ could help to determine these 
form factors. 

In nuclei, the description of pion production requires a realistic model to account for the final state interaction 
(FSI) of the pion. This is usually implemented in MonteCarlo codes. However, that is not enough. It has been shown 
in several works that also the production mechanisms are modified in the medium. Furthermore, in some particular 
cases, like the coherent pion production, a quantum treatment of FSI is neces sary . In any case, there is clear progress 
in our understanding of this reaction and some recent [25| and coming data [26| | may already put strong constraints 
on the A form factors [27| . 

Above these energies, but still below the DIS region, new inelastic channels are open and several baryonic resonances 
beyond the A(1232) can be excited El [28j] . Recently, there has been an important progress in the determination of 
their vector form factors with the advent of high quality electromagnetic data [25^ |301.[3~l11 . Our knowledge of the axial 
form factors is, in general, poorer due to the scarcity of experimental information [1S| . 

The first of these resonances is the A''* (1440), for which the weak excitation can allow to study the axial sector. 
Although its main decay channel is to A7r, its contribution to the one pion production cross section has been found to 
be negligible [3] because of the strong dominance of A mechanisms. However, the situation can be different for the 
production of two pions. This channel starts at invariant masses of the hadronic sector just below the A. However, 
the A does not couple to two pions in s-wave and thus it is not very relevant at these energies, where only slow pions 
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are produced. On the other hand, the Roper resonance TV* (1440) has a sizable decay into a scalar pion pair and it is 
very wide so that its contribution could be large. Indeed, Roper excitation mechanisms are known to play a major 
role in other two pion production reactions close to threshold, like ttN — > ttttN [H, [H, [13] or N N — > nnNN [H, [3(| . 

In the weak interaction sector, there exist very few attempts to study the two pion production induced by neutrinos 
and antineutrinos. The older experiments done at CERN [37|, |3|| |39|, |4Cj , in the regime of high energy have studied 
two and three pion production to investigate the diffractive production of meson resonances like p and rj. The 
later experiments done at ANL [HI , |42T | and BNL [43[ at lower energies have investigated the two pion production 
processes, specially in the threshold region, in order to test the predictions of the chiral symmetr y o f the strong 
interaction Lagrangian. Such studies were theoretically proposed by Biswas et al. [44j and Adjei et al. [451 146| . Biswas 
ct al. used PCAC and current algebra methods to calculate the threshold production of two pions. On the other hand, 
the work of Adjei et al. made specific predictions for the threshold production of two pions in a restricted kinematic 
region using an effective Lagrangian incorporating chiral symmetry and its breaking, governed by a free parameter 
(£). Imposing these kinematical restrictions, the experimental data of ANL and BNL were analyzed and compared 
with Adjei et al. results. However, the model did not include any resonance production and the contribution of the 
A* (1440) should be taken into account. Furthermore we use an expansion of the chiral Lagrangian that includes terms 
up to 0(1/ f 3 ), while Adjei et al. kept only terms up to 0(1/ f 2 ). Now that the pattern of chiral symmetry breaking 
is well known and the background terms contribution to the threshold production of two pions is fully determined, 
the process could be used to study the electroweak transition form factors of the Roper resonance. 

In this paper, we will study the viN — ► 1~-k~kN channel close to threshold. Apart from Roper resonance contribution, 
many other background mechanisms, that only involve nucleons and pions, appear and are described by an effective 
Lagrangian. We will use the lowest order chiral perturbation theory Lagrangian to derive the needed axial and vector 
currents. 

In Sec. [Til we present the formalism and the Lagrangians used in our model. We also give the expressions of 
the Roper form factors. In Sec. IIII1 we present our results and compare them with the available data. Finally, the 
appendix gives the detailed formulas of the contributions of the background mechanisms for all channels. 



II. MODEL FOR v INDUCED TWO PION PRODUCTION 



A. Kinematics 



We will focus on the neutrino-pion production reaction off the nucleon driven by charged currents, 

v x (k) +N(p)^r(k')+ N(p') + n(k %1 ) + n(K 2 ) (1) 

though the generalization of the obtained expressions to antineutrino induced reactions is straightforward. The unpo- 
larized differential cross section, with respect to the outgoing lepton kinematical variables, is given in the Laboratory 
(LAB) frame by 

dQ(k')dE'-^ |£|^^cc 2 J (2) 

with k and k! the LAB lepton momenta, E' — (k 12 + m 2 ) 1 / 2 and mi the energy and the mass of the outgoing lepton, 
G = 1.1664 x 10~ n MeV -2 , the Fermi constant. We take e i 2 3 = +1 and the metric g^ v = (+,-,—,—), thus the 
leptonic tensor is given by: 

L^cr — (L s )^ l(T zh i(L a )p iCT — k^k a + k a k^ — g^crk ■ k zh is^aapk a kP (3) 

where the +(— ) sign corresponds to neutrino (antineutrino) induced processes. 
The hadronic tensor reads as 

Wg^ = 

£ / ^^ind^2^ (4) 

spins v ' JV\/ iv/ 

with M the nucleon mass, q = k — k' and E' N the energy of the outgoing nucleon. The bar over the sum of initial 
and final spins, denotes the average on the initial ones. As for the baryon states, they are normalized so that 
(p\p'} = (2 7r ) 3 S 3 (p — p')po/m. By construction, the hadronic tensor can be split in 

(W^ 2v ) = (W^) s + i(W^) a (5) 

with (Wqq 27T ) s and (W^ 27T ) a real symmetric and antisymmetric parts, respectively. 
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B. Lagrangians for the non resonant terms 



For the derivation of the hadronic tensor we use the effective Lagrangian of the SU(2) non-linear a model. This 
model was used previously in [2l| for the description of the non-resonant contributions to one pion weak production 
processes off nucleon. We refer the reader to that paper for details. Up to 0(1/ f 3 ), this SU(2) chiral Lagrangian 
reads 



£ = *[$- A/]* + \djd^ - \m 2 J 2 + £f nt 



(6) 



C 7 



< !A y^ 75 l(dJ)V - -L§ 7A1 f (0 x * - i ($%$B*$- ($d^)($d»4>) 



f; 



9A 

'6/1 



$ 2 ^dJ-($dJ)^-$ 



l TT (±2\2 



24/2 



(7) 



where VP = 



P 



is the nucleon field, <fi is the isovector pion field, f are the Pauli matrices and f„ = 93 MeV is the 



pion decay constant. The vector and axial currents generated from the Lagrangian in Eq. (J6j) are given by [2l| 



y fl ^ v » ' 



*-^($*d^)+0(±) (8) 



■I 7T 



•0*V) -0 2 ^(/) 



*+o( 75 ) 



(9) 

and determine the weak transition vertex where the W— boson is absorbed. These currents are, up to a factor, the 

Thus, A£ 
vector and 
contribute 



hadronic realization of the electroweak quark current j^ c , for a system of interacting pions and nucleons 



and VJf account for the W— decay into one and two pions, respectively, while V/ 1 and A 1 / provide the WNN ■ 
axial vector couplings. Besides, A% and lead to contact WNNtt vertices and finally A^ and V// either 
to processes with more than one pion in the final state or provide loop corrections. 

The overall normalization can be obtained, for instance, by relating the currents of Eq. © and Eq. @ with the 
phenomenological vector and axial nucleon currents in the (-/V , 7r|j£ 1 c+ (0)|./V) matrix element 

fop' =p + q\j« + (0)\n;p)=co S & c u(p ')(V%(q) - A%(q))u(p) = A a (10) 

where the u's are Dirac spinors for the neutron and proton, normalized such that uu = 1, and vector and axial nucleon 
currents are given by 



F 2 V (q 2 ) av 



A%(q)=G A (q')x 7 Q 75 + 



jq 75 



(ii) 



We find 1 that -y/2 cos8c ([V^+i — provides the Vy + -absorption vertex, with the appropriate normalization. 

The magnetic part in Eq. pip is not provided by the non-linear sigma model, which assumes structureless nucleons. 
We will improve on that by including the q 2 dependence induced by the form factors in Eq. (|11[) and adding the 
magnetic contribution, F% term, to the vector part of the W + N — > N amplitude. 



1 The +1 spherical component of a vector A is defined as A+i = — (A x + iA y ) /y/2. 
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1) m) n) 

FIG. 1: Nucleon pole, pion pole and contact terms contributing to 2n production. 
For the nucleon vector form factors (FF) we use the parameterization of Ref. [56| 



N (j e j- ku m N u- M e rip b i/ ~*M r-i , \ ~\~E I ± \ fio\ 



with k = -q 2 /4M 2 , M D = 0.843 GeV, (i p = 2.792847, /u n = -1.913043 and A„ = 5.6. 

F?(q 2 ) = \ {Fl{q 2 ) - F?{q 2 )) , ^(q 2 ) = \ (/^ 2 V) - H n F?(q 2 )) . (13) 
The axial form factor is given by (57j 

GA(q 2 ) = j, g /L. 2 , 2 , Sa = 1.26, M A = 1.05 GeV (14) 

(1 - q 2 /M 2 A ) 2 

Using this current we obtain the sixteen Feynman diagrams, depicted in Fig.[l] constructed out of the W + N — > N, 
W + N — > Nir, W + N -> 7V7T7T, and the contact W + ir — > 7r weak transition vertices (Eqs. (JEHU) and the irNN, irirNN, 
wirirNN, tttttttt couplings (Eqs. (J6HZ])). Since we have included a q 2 dependence (F± (q 2 )) on the Dirac part of the 
vector WNN vertex and to preserve vector current conservation, we include the same form factor in the Vr b d weak 
operators. 

The expressions for the matrix elements of these terms, for all possible channels, can be found in the appendix. 
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FIG. 2: Direct (left) and crossed (right) Roper excitation contributions to 2tt production. 



C. Contribution of the iV*(1440) resonance 

The Roper iV*(1440) P\\ is the lowest lying baryon resonance with an s-wave isoscalar two pion decay. This suggests 
the possible relevance of Roper excitation mechanisms in two pion production processes close to threshold. Indeed, 
its importance has been clearly established in the nN — > ttttN [32|, [33| and the NN — » Ntttt [35L [36[ reactions where 
it plays a dominant role for certain isospin channels. However, that is not the case for electromagnetically induced 
reactions, due to the relatively weak coupling of the Roper to the photons. See, for instance, Ref. [471 ] . 

We include in our study the two mechanisms depicted in FigJ^l which account for the Roper production and its 
decay into a nucleon and two pions in a s-wave isoscalar state. However, we do not include the contribution of 
mechanisms in which the Roper decays into 7rA(1232) states, because the two produced pions are in p-wave and thus 
are not expected to be relevant close to threshold. 



1. Roper interaction Lagrangians and currents 



The Lagrangian for the s-wave N* — > Ntttt decay can be written as 



neglecting terms of order 0(p 2 /M* 2 ) [48j]. Here, ipN* is the Roper field. In Ref. [35| the best agreement with 
NN -> ttttNN and the ttN -> ttttN data was obtained using c\ = -7.27 GeV -1 , c* 2 = GeV -1 . This result was 
obtained assuming a branching ratio of 7.5% for the N(tttt)jZ^ decay mode and a total decay width of the N*, 
r tot = 350 MeV gg|. 

The other required ingredient is the coupling of the Roper to the charged weak current, that is the vertex W + n — > 
7V*+(1440). The matrix elements can be written as 

(N*+;p* = p + q\j? c+ (0)\n;p) = cos0 c u,(fc) J^u{p) , (16) 

where it* is the Roper spinor and 

Jcc* = 2 — W i - H + l ~ & Qv - Gat 75 3 ?75 <* qvls (17) 

I* H MM 

is the most general form compatible with conservation of the vector current and Dirac equation for the nucleon and 

Roper Dirac spinors. A factor fi = M + M*, with M* = 1440 MeV the mass of the Roper resonance, is introduced 

in order to make the vector form factors dimensionless. The Gt term, and unlike the elastic nucleon case where it is 

zero due to G-parity invariance, does not need to vanish in the present case because the nucleon and the Roper do 

not belong to the same isospin multiplet. Nevertheless, most analyses neglect its contribution and we shall do so here. 

The N*Nir coupling is described by the pseudovector Lagrangian 

Cn*n* = — *at-7 m 75t • + h.c. . (18) 

The decay width for this process is given by 

T N ,^ N = ^- (-L) ^|g cm | 3 , (19) 
27T \ m,r / W 

where W is the N* invariant mass and \q cm \ is the momentum of the outgoing pion in the outgoing ttN center of mass 
frame. Taking for this decay channel a branching ratio of 65% and the total width of the Roper T — 350 MeV [49| . 
we get / = 0.48. 
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2. N* (1440) form factors 



Assuming the pseudoscalar coupling Gp is dominated by the pion pole contribution, and imposing partial conser- 
vation of the axial current (PCAC) hypothesis we can relate Gp with the axial coupling Ga 



G P (q 2 



m 2 — q 2 



G A (q 2 ) ■ 



(20) 



Furthermore we can relate Ga with the N*Nir coupling constant at q 2 = using the non-diagonal Goldberger-Treiman 
relation 



G A (0)= 2^^=0.63 



(21) 



The q 2 dependence of Ga(<Z 2 ) is not constrained by theory so we shall assume for it a dipole form of the type 



G A (q 2 ) = 



G A (0) 



(i- q 2 /M 2 A j 2 > 



(22) 



with an axial mass Ma* = 1 GeV. 

The vector- isovector form factors FY* / '/i 2 and FY* / h can be related to the isovector part of the electromagnetic 
(EM) form factors, FY* — Ff* — FJ 1 *, that can be determined from photo- and electroproduction experiments. The 
relevant experimental information about these EM form factors is usually given in terms of helicity amplitudes for 
the EM current j£ m , defined 2 as [5l[ 



pol 



Sfo = J^-j^=(N* t | je.nu(0)|JVt>£. 



kit ^J-q 2 



(23) 
(24) 



pol 



where N stands for proton or neutron, a = 1/137, q is the momentum of the virtual photon, kp = (W 2 — M 2 )/2W, 
with W the energy of the Roper in its center of mass, and the polarization vectors are given by 



e ± = -±=(0,Tl,-*,0), 
and for a photon of momentum q moving along the positive z-axis 

' KM, 0,0,?°). 



(25) 



(26) 



The factor £ in Eqs. (|23|) and (|24[) is given by the relative sign between the NNir and N*Ntt couplings |51[ which we 
have taken to be positive (we will discuss this point later). 
Finally, the EM 7 TV — > N* current is written as 



(N*;& =p + q\£ m X0)\N;p) = 
Using Eqs. (p?4"l) - (f2"7| we obtain the following relations 



</ 2 = l?W 



u{p) . 



r t?n* 



q 2 F< 



N* 



H W + M p 2 



Ft 



N* 



?N* 



1 



\i W + M 



(27) 

(28) 
(29) 



Please note that the MAID group analysis [50| define with the opposite sign, which we take into account when comparing to their 
data. 



7 



with 

2 = / 8na(W + M)W 2 
9[q ' y M(W-M)((W + M) 2 -q 2 )' [ ' 

Inverting Eqs. (]28[) and (|29|) we can obtain the proton FF Ff* as a function of the experimental helicities. Unfortu- 
nately, there is no such information on the 777, — > iV* transition so we need some theoretical assumptions to relate 
the isovector helicity amplitudes with the EM ones. Following the quark models predictions of [53| and [54[ we shall 
assume A^ 2 = —2/3A^ 2 and 5"y 2 — 0. For a review of different models see Ref. [55[. Using these relations we can 

write the neutron FF as a function of the proton ones and express the vector FF in terms of only Ff and F 2 * as 

v , _ FTm + Wf - 5g 2 /3) + 2/3Ff(M + gQg 

v , _ Ff(5(M + HQ 2 - 3g 2 ) M - 2JfV(M + TU) 

2 ~ 3((M + IU) 2 -g> • 1 1 

We have fitted the proton-Roper EM transition form factors to the experimental results for helicity amplitudes 
given in [30l l3l| using a parameterization inspired by Lalakulich et al. 18] 



1 - <?/X x M* 

rP / / « 2 



jtvh , ^T,^ , (33) 



Ff(g2) = ^_ x _ X2ln ! _ ^_ ) ) (34) 



D v V V IGeV 

where = (l — q 2 /My) 2 with My = 0.84 GeV. We have fitted the dimensionless parameters g\ , g%, Xi and X2 to 
the available experimental data, and found the following best fit (labeled as FF 1 in the results): 

5 [ = -5.7±0.9, gl = -0.64 ±0.04, X 1 = 1.4±0.5, X 2 = 2.47 ± 0.12 , (35) 

with a x 2 /d.o.f = 5.2. Our best fit parameters g\ and X\ differ appreciably from those quoted in Ref. [3], in 
particular g\ comes out with opposite sign. This is because the authors of this latter reference did not consider the 
existing extra minus sign among their definition of which agrees with that of Eq. I|24p and used in this work, and 
the definition used in the MAID work. Indeed, if we do not consider this minus sign, we find values of the parameters 
in good agreement with those quoted in with minor differences due to the use of different data sets. 

The Roper EM data have large error bars and it is possible to accommodate quite different functional forms and 
values for these FF. Thus, we shall compare other different models for the vector form factors. Firstly, we consider 
the constituent quark model with gluon, pion and tr-meson exchange potentials as residual interactions of Meyer et 
al. [E^l in set FF 2. In this model the electromagnetic current included, in addition to the one-body current, two-body 
exchange currents associated with the quark-quark potentials. Here no assumption about the relation between proton 
and nucleon form factor was assumed. We will also consider the recent parameterizations 3 of Lalakulich et al. [IH in 
the set labeled FF 3. Note that this work employed a pseudoscalar form for the irNN coupling used in deriving the 
Goldberger-Treiman relation, instead of the pseudovector, Eq. (|T5]) . We shall finally use the predictions of the recent 
MAID analysis [H in the set labeled FF 4. 

From all this information we can now obtain definite expressions for the currents of the direct (N*d) and crossed 
(N*c) diagrams shown in Fig. [5] 

A%, d = 2g*u(p')S*(p + q)J^(qHP) (36) 
A^ c = 2g*u(p')j^(q)S*(p' - q)u(p) , (37) 

with g* — cl(m v / fjr) 2 . In the crossed diagram it appears J^ c *, that corresponds to the crossed vertex W + N* — * N 
and it is given by Jc„ = 7°(J^c*) t 7°- We have here introduced the Roper propagator 

S * iP * ] = P l-Ml+!{lt+*W)T tot {W)l2 ■ (38) 



3 Notice the opposite sign convention for Fy in 18] 
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The total W-dependent decay width rtot(W^) includes all the possible decay channels, namely Nir, given by Eq. (|19p . 
7rA and N(tttt)^q. For these two cases we take 

T w a(W) = 350 x 0.275K(WO| 3 /K(M*)| 3 McV (39) 

and 

T^n = 350 x 0.075PhSp(W)/PhSp(NQ MeV , (40) 

where PhSp(W) is the phase space for the three body 7T7riV decay, and p^ is the pion momentum in the Roper rest 
frame. 

The isospin structure of Eq. (|15p determines that the only allowed channels for the two pion s-wave decay of the 
Roper resonance are tt + tt~ and 7r°7r°. The above expressions are given for the tt + it~ channel, therefore we must 
include an additional 1/2 symmetry factor for the 7r°7r° channel. 

3. Relative sign between the Roper and non-resonant contributions 

Here we give some details on the used scheme to set up the relative signs between the different contributions: 

• All relative phases in the non-resonant terms are completely fixed by the non-linear sigma model Lagrangian of 
Eqs. © and ([7|), the currents deduced from it and shown in Eqs. © and ©, and the phenomenological WNN 
vertex. 

• We have assumed the sign of the N*Nn coupling to be the same as that of the NNtt coupling (positive). This 
choice fixes the global phases in Eqs. ([23| and (|24|) [5lj ] and thus the phase of the vector part of the WNN* 
current. Besides the non-diagonal Goldberger-Treiman relation (Eq. ([2Tj) ) fixes the axial part of the WNN* 
current. 

• Furthermore, once the relative sign of the NNir and N*Nir couplings has been set, the study of the reactions 
NN — > NNtttt and irN — > Nirn in Ref. [3|| fixes the values for the c\ and c\ N* Nitty couplings. 

Taking the opposite sign for the N*Nir coupling, does not affect the results. This is because the signs of both, 
the WNN* current (see previous discussion), and the c* and c*, couplings 4 would change and therefore the whole 
resonant contribution would not be affected at all. 



III. RESULTS 

In this section, we present the results for all the two pion states produced in the neutrino induced reactions on 
nucleon targets, i.e. vp — ► u~pTr + ir°, vp — > n~nir + ir + , vn — > n~pir + ir~ , vn — > n~nir + ir° and vn — > u~pir°ir . 

Our model includes all relevant terms for neutrino induced two pion production close to threshold. In addition to 
the contribution of nucleon and pion pole and contact terms described by the chiral Lagrangian the contribution of 
the iV*(1440) resonance coupling to two s-wave pions is included. The A(1232) and other higher resonances are not 
considered as their contributions would vanish at threshold. 

In Fig. [31 we present the results for the cross section for the process vn — ► u~pTr + ir~. We show separately 
the contribution of the background terms coming from the nucleon pole, pion pole and contact terms as well as 
the contribution of the Roper resonance as calculated by using the various form factors described in section [Til 
The interference between background and the Roper contribution is not shown. We see that the background terms 
dominate the cross section for neutrino energies E v > 0.7 GeV. At lower energies the contribution from the Roper could 
be larger or smaller than the background depending upon the vector form factors used for the W + NN* transition. 
The parameterization determined by the recent MAID analysis 50] gives the largest Roper contribution to the cross 
section. The differences in the predictions for the cross sections using the various parameterizations could reach a 
factor two. The Roper contribution is specially sensitive to F)f*{q 2 ) which is negative in contrast to the positive value 
which one gets in the case of the nucleon. This has interesting implications for the vector-axial vector interference 



4 These couplings would swap sign when the sign of the N*Nn coupling is changed, since they are fitted to the NN — » NNtttt and 
ttN — > Ntttt data. 
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0.6 0.7 0.8 0.9 

EfGeV) 



FIG. 3: Cross section for the vn — > /i pn + n reaction as a function of the neutrino energy. The interference between 
background and the N* contribution is not shown. See text for details. 



contribution for the cases of v and v excitation of the TV* (1440). The comparison of data on v and v induced two 
pion production in this channel, if available in the future, could be used to study the WNN* transition and better 
constrain the corresponding form factors. 

In Fig. [U we compare our results for the same process with experiment. The largest results are obtained with the 
set 4, and the smallest from the set 3 of nucleon- Roper transition form factors. The other two sets give cross sections, 
not shown in the figure, lying between these two limits. Our results with the full model are larger than those obtained 
using only the background terms as done by Adjei et al. [Ifl, however the results are still lower than the central value 
of the experiment. Experimental data with better statistics is highly desirable in this channel to make a decisive 
comparision with the predictions of our model. 

We expect our model to have a limited region of applicability from threshold to an invariant mass of the ttttN 
system below 1.4GeV. For muon neutrinos, that implies a LAB energy under 750 MeV, beyond which some additional 
contributions from the A(1232) and other higher resonances will become relevant. In order to make a meaningful 
comparison with the theoretical calculations for threshold two pion production, Adjei et al. [451 ] sug gested a kinematical 
cut on phase space which was implemented in the experimental analysis made by Day et al. [421 ] and Kitagaki et al. 
[43l | . These cuts were defined as 

ql < ((l + n/2)m 7r ) 2 , (41) 
P ■ q-K < (M + (1 + v)m^ f - M 2 - m\ (42) 
p' ■ < (M + (1 + V )m^) 2 - M 2 - m \ , (43) 

with qv = (k ni + fc 7r2 )/2. Different choices for 77, specifically 77 = 1/4, 2/4, 3/4, were proposed. As explained in 
Ref. [451 ] . the first inequality keeps the individual pion momenta close to the average pion momentum, while the last 
two restrict the phase space to regions of small invariant mass of the three hadrons. 

Using these kinematic restrictions on the phase space with rj = 3/4, as used by Refs. [42|, [43j], we present the 
results for the cross section for the vn —f (i~pTr + Tt~ channel in Fig. O and compare with their data. We show our 
results with only background terms and with the full model evaluated using the set 1 of nucleon-Roper transition 
form factors. Other sets give a similar result in this case. We find that, even in this kinematic region, the theoretical 
results including the resonance contribution are lower than the experimental. 
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In Fig. [SI we present the results, with the same cuts, for the total cross section for the channel vp — > fi~mr + Tr + 
and compare with the data of Ref. (42J. For this channel there is no contribution from the iV*(1440) resonance. Our 
results are in agreement with those of Adjei et al. [45j ] and are consistent, within errors, with the experiment in 
the higher energy region. However, we underestimate the first point (E v = 1.25 GeV) by more than one order of 
magnitude. This d isag reement may be ascribed to the low statistics of the experiment after the kinematical cuts have 
been implemented [42| , to inadvertant inclusion of some P wave pions while implementing the kinematical cuts in the 
experimental analysis but also to possible additional reaction mechanisms not included in our model. Note that the 
kinematics of the plot explore the region just up to 1.4 GeV of final hadron state invariant mass. At this energy we 
do not expect our threshold production model to account for all the mechanisms of pion production, many of which 
would be in p-wave and therefore beyond the scope of our work. This mechanisms would include A resonance terms 
similar to those studied in [47| in the context of photo-induced reactions. In any case, the disagreement at E = 1.25 
GeV underscores the need for better experimental data in this channel around neutrino energy E — 1 GeV. 

Finally, in Fig. [2 we show our results for the threshold production of two pions in the other channels. Again we 
only show our results with only background terms and with the full model evaluated using the set 1 of nucleon-Roper 
transition form factors. The use of different sets gives similar results for the full model calculation. There is no data 
in the low energy region; the limited data available in these channels at higher energies have not been analyzed in 
the kinematical region of the threshold production of two pions and thus a comparison with our present results is not 
possible. 

To summarize, we have studied the threshold production of two pions induced by neutrinos on nucleon targets using 
a chiral Lagrangian to calculate the nucleon pole, pion pole and contact terms. The contribution from the excitation 
of the Roper resonance, A^*(1440) has been included. The vector transition form factors are determined from the 
available data on the helicity amplitudes of its electromagnetic excitation. The axial form factors are obtained using 
PCAC and the experimental data on N*Nw decay. The q 2 dependence of the axial form factor has been assumed 
to be of the dipole form. It is found that the Roper resonance contributes strongly to the vn — > /i~p7r + 7r~ and 
vn — > fj,~pir°ir° channels. Its contribution is as important as the contribution of the nucleon, pion pole and contact 
terms, up to energies of E v ~ 700 MeV. The resonance contribution is sensitive to the form factors used for the WNN* 
transition, and the experimental data, if available with better statistics for the vn — > n~pir + ir~ and vn — > /i _ p7r°7r° 
channels, could be used for the determination of these form factors. For other channels, as well as for these, the 
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E v (GeV) 

FIG. 7: Cross sections as a function of the neutrino energy. All calculations correspond to the full model with the FF1 set of 
nucleon-Roper transition form factors. 

theoretical results underestimate the experimental results. Furthermore, physical channels involving A production 
could become relevant at beam energies higher than 750 MeV, where the threshold approximation we have assumed in 
our model would no longer be valid. Nevertheless, the low statistical significance of data does not provide a conclusive 
test of our model. Availability of improved data, specially from future experiments on neutrino induced two pion 
production in the region of energies E v < 1 GeV will be very useful in the study of the electroweak properties of the 
TV* (1440) resonance. 
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APPENDIX A: HADRONIC CURRENT FOR THE DIFFERENT CHARGE CHANNELS 

The contributions are labelled following Fig. [TJ We give explicit expressions for all channels except for W + n — > 
pir°ir° which can be obtained via the following isospin relation 

A{W + n -> pir°ir ) = A{W + n -> pn + n~) - (A{W + p -> pn+ir ) - A{W + n -> ?i7r+7r )) . (Al) 

V2 

In all these expresions we have used for the nucleon propagator 



13 



and — are denned as in Eq. (fTTj) . 

1. W+pp -^Pp'7r^7rg 3 



At = cos 9 C u(p') (3fe + 2M ) 75 u (j5) (A4) 



.41' = ^cosfl c ? f— % ^ 24/77, i 7',„„!/7 ) (A ,, 



A» = cos C {q _l 2)2 q _ m2 2FY W)h »(0 (AT) 



<' = -tS cos ^ "(P') ( V " - j4 ") 5 (p' - ?) (ft - fc) (A8) 

J 7T 



^ = cos C (fe 7 5S(p' + fe 2 ) 2*JWtS - J 5 S(p' + k 2 )^ ) u(p) (A9) 

^ = COS0C {W)^g A ^75 (S(p - k 2 )fc lb + 2S(p - h)lt l75 ) u(p) 

-u(p')r (- S( P - k 2 )hl5 + 2S(p- feigns) u(p)} (A10) 

^ = --^-r cos0 c 2 u(p') fe7 5 S(p' + fca)« + (All) 

At = --^-cos6 c g " u(p') (d + h)S{p-k 2 )h lb ~2{i^h)S{p-k 1 )h lb ) u(p) (A12) 

a 2 — 9k^ 

A J = dfe COS0C 2F i ( n \v 2 W)h^S{p' + k 2 ){i - ^) 75 u (j5) (A13) 

n 2 ( — 2fr M 

At = cos 9 C 2F? { (g,^^. 1 ^ - W75S(P - fc 2 )^ 275 «® 

- 2&2 1 



(q - k 2 ) 2 - m 2 



u{p'){i - fch 5 S( P - u(jO (A14) 



^ = --4- cos^ c u(p')h^S( P ' + k 2 ) (V" - A**) S(p - fci)^ 75 u(p) (A15) 

= inhi cos9c ( V " - A ^ S( -P - kl ~ fc2 ) (^75^(P - fci)&76 - #175% - fc 2 )^ 275 ) (A16) 

2V2/ 2 
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2- W+p p -> n p /7T+7r+ 



A% = — cos C u(p') r - 9 a 7^75) u(p) (A17) 

/•7T 

A = £72 cos ^ 2 ^ 2 ( 3 ^i + 3 ^ + 4M) 75 (A18) 

"Jot 9 TO ir 



3A g M f (4gfci - 2 g fc 2 + 2/dfc 2 - g 2 ) (4gk 2 - 2gh + 2k x k 2 - q 2 ) \ , , .., _ 

^ - 2^ J 2 cosO c { {q % q l m f &)hum + (g _ 2 t)^^ ^"Cg) > (A2D 



2/2 

(^tsS(p' + fci) + ^2 J 5 S( P ' + k 2 )) -futf)} (A22) 



A f = -^k C0S ^ { (^75^(P' + fci) + 7 5 S(p' + **)) 2f\ y 3^75 u(0) 



Z Jtt 

-u{p')Y (S( P - fci)^i 75 + S(p - fc 2 )^275) u(p)} (A23) 
-< = -772 cos ^ 2^ 2 {h-nStf + fci) « + #2) + &7 5 S (p' + k 2 ) (i + ft)) u (0) (A24) 

4 /7r 9 TO 7T 



A- = TH cos ^ 2^ 2 W) (« + h) S(P - fc 2 )^275 + d + h)S(p - fci)^i75) «(£) (A25) 



„2 / n n _ out 1 



+ (g ! fc2 )2_ 2 m 2 W)h^S{ P ' + k,){i fe)75«(^J (A26) 



2 cos W - \ 1 2 u(p')(rf - hhsS(p- k 2 )f 2l5 u{p) 



fl^ — 2k^ 

+ {q -k 2 ) 2 -ml fchsSip-h^uip) ) ; A27) 



cosfl c (^i7 5 5(p' + fci) (V" - A M ) S(p - fc 2 )^ 75 + fc^S{p' + k 2 ) {V» - A 11 ) S(p - fci)^ 75 ) u(p) 

(A28) 
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3. W+ 



A£ = ~ cos 6 C u(p') {2FY r - g A 7^75) u(p) (A29) 

1 Jtt 



A b=~7rk cos c 1^ 2 ( 3 #i + 2M ) 75 u(p) (A30) 

6/ 2 - mi 

9n a — 3£- M 
f (4gfc 2 - 2gfcj + 2fcifc 2 - g 2 ) , 3m 2 



cos0c ?^ 1 - ( ; - fc ; = k2 y - mi J + {q - fcl - - ro » ^ ^ ^ 



1 2£ M — 

A " = -j2 COs9 c (q-k^-ml 2F ^(p)h u(p) (A33) 



A a' = 772 cos °c W) (h - h) S( P + g) (V" - A^) u{p) (A34) 

A% = ~ cos6 c u(p>) (V» - A^) S(p' - q) (h - fr) u(p) (A35) 

A f = -§72 cos °c W) {h^S(p' + k 2 )2F^g A ^ l5 - fi 2l5 S(p' + k 2 )r) u(p) (A36) 

A 9=-^72 cos 9 * (2F 1 y 5A7 M 75^(p ~ ^2)^275 - -fS{p - fc 2 )^ 7 5 ) «($ (A37) 



A h = jk cos0 c «(p')#27 5 S(p' + fc 2 )« + (A38) 



A t = -772 cosd - 2!" 2 + - fc 2)^275 u($ (A39) 



^ = ~2% COS0C 2F ? ( q l kl y-ml ^P'W^iP' + M« - ^i)75 u{p) (A40) 

a 2 — 9k^ 

A l = ^k C0S ^ 2F i ( -TY2 1 2 "GOM - ftb*S(p - fc 2 )^ 75 «® (A41) 

A ? = -^H cos ^ W)h^S(p' + k 2 )hlsS{p' + ki + k 2 ) (V - A* 1 ) u(p) (A42) 

A n = -^ cos6 c u(p') (V» - A^) S( P -h- k 2 )^ 5 S(p - k 2 )hl5 u{p) (A43) 
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V^Tlp -> Tip' 7T^ 7Tfc 2 



K = -^Jl cos 0c (2^i V 7 M - 9 A 7 M 7 5 ) (A44) 



= "rfe cos 0C 2" 2 «(P') ( 3 + 2M ) 75 u(p) (A45) 



^ - -^-^ fa,^.^ 1 ,^ 2M^«© (A46) 

< = 7=^ costf c 4^-2^ + 2^^-^ 2M _ 

d 3V2/2 g 2 - m 2 . (g - fa - fa) 2 -ml 

1 2fc M — o M 

^ = " V2p C ° S9c {q-kf -ml 2F ^(p)h u{p) (A48) 



A *> = -li cos0 - W) (fa ~ fa) S(P + q) (V" - A") u(p) (A49) 



A) = -^^cos9 c { u(p') (fa l5 S(p' + fa) ~ 2fa l5 S(p' + fa)) 2F^g A ^ l5 u(p) 

-u(p') {faj 5 S(p' + fa) - 2fa l5 S( P ' + fa)) ^ u(p)} (A50) 

A ° = ~27Tp COS9c ( 2F i9A^7 5 S(p - fa)fa lb - ^S(p - k 2 )fa l5 ) u(p) (A51) 
At L = --2=rr cosg e 2 ^ 2 u(p') (faj 5 S(p' + fa)(i + fa) - 2fa l5 S{p' + fa)(4 + fa)) u(p) (A52) 

At = --7=tt cos^ c 2 + fa)S{p - fa)fal 5 u(p) (A53) 

4V2/2 q- - mi 

n 2 / — Ik^ 

Aj = cos0 c 2F x y ^ {q l kl) 2l m 2 u(p')faToS(p' + k 2 )(4 - fa) 75 u(p) 

-2 , \ , 2 2 2 V(?)h>ftSW + fa)(4 - fc)j 5 u(p) (A54) 
{q - fa) 2 - mi ) 

n 2 — 9k^ 

At = ^- 2 cos6 c 2J^_L_^_ - ^) 75 S(p- fa)fa^u($) (A55) 

^ = cos c «(p ') (fe7s5(p' + fc 2 )#i75 - ^i75^(p' + fa)fa l5 ) S(p' + fa+ fa) (V - ^) (A56) 

A£ = cos^ c Vtf'Wx'IsSijf + fa) (V - A") S(p - fc 2 )fe 75 u(p) (A57) 
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